Electron-electron interactions in the conductivity of graphene 
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The effect of electron-electron interaction on the low-temperature conductivity of graphene is 
investigated experimentally. Unlike in other two-dimensional systems, the electron-electron inter- 
action correction in graphene is sensitive to the details of disorder. A new temperature regime 
of the interaction correction is observed where quantum interference is suppressed by intra-valley 
scattering. We determine the value of the interaction parameter, Fq — —0.1, and show that its 
small value is due to the chiral nature of interacting electrons. 
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The low-temperature behaviour of the resistance of 
electron systems is determined by quantum effects. Two 
distinct phenomena are responsible for this: quantum in- 
terference of electron waves scattered by impurities (weak 
localisation) WL [1| , and electron-electron interaction in 
the presence of disorder, EEI Q. The WL correction is 
used to study electron dephasing, while the EEI correc- 
tion, which is not sensitive to dephasing, has been used 
to probe the dynamics of interacting electrons, e.g. [H-dj]. 
In graphene, the charge carriers are chiral and located in 
two valleys. As a result, WL is sensitive not only to in- 
elastic (phase breaking) scattering ,but also to a number 
of elastic scattering mechanisms [3, |lfj|. For this rea- 
son, the WL correction to the conductivity can be either 
negative or positive, depending on the experimental con- 
ditions 
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So far the effects of EEI on the low-temperature con- 
ductivity of graphene have not been studied experimen- 
tally, and only the high-temperature ballistic regime was 
analysed theoretically (l3j ]. It was predicted that at 
ksTTp > 1, where t p is the momentum relaxation time, 
the EEI correction is determined by coherent backscat- 
tering on a single impurity, which in graphene is sup- 
pressed due to the chirality of charge carriers [l4| . As a 
result, the EEI correction can only occur due to scatter- 
ing on atomically sharp defects and is expected to have a 
universal form which is independent of the details of the 
electron-electron interaction. In addition, the ballistic 
regime in graphene can only be realised at relatively high 
temperatures where the EEI effect has to be separated 
from strong effects of electron-phonon scattering. Thus 
the ballistic regime of EEI is not expected to be promis- 
ing for the study of electron interactions in graphene. 

In the diffusive regime, kgTT p < 1, interacting elec- 
trons in graphene scatter on multiple impurities, so that 
backscattering is less important. Hence the EEI correc- 
tion will contain information about the details of inter- 
action. In this work we study the EEI effect on the con- 
ductivity of graphene in this regime. To separate the cor- 
rections due to EEI and WL, we combine measurements 
of the temperature dependence of the conductivity with 
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FIG. 1: (Color online). The gate voltage dependence of the 
resistivity for samples SI, S2 and S3. The bars (shown for 
S2) indicate three studied regions. Insets: graphene band 
structure with two valleys, and a schematic illustration of the 
regimes of interaction between two electrons that are scat- 
tered by impurities. 



studies of magnetoresistance (MR). We show that our 
results are described by a logarithmic correction to the 
conductivity [2j: 



5cr (T) = -A(Fg) 
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Here the coefficient A(Fq) is determined by the strength 
of interaction and the symmetry of electron states. We 
show that in graphene there is a new regime of EEI and 
find the value of the interaction parameter Fq . 

Three samples with Hall-bar geometry (SI, S2 and S3) 
are fabricated by mechanical exfoliation of graphite on 
Si/Si02 substrates (ro| . Sample parameters are shown 
in Table 1. Quantum Hall effect measurements were per- 
formed to verify that the samples are monolayers [15]. 
Figure Q] shows the resistivity p as a function of the gate 
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TABLE I: Electron mobility fi (in cm 2 V _1 s _1 ) for three re- 
gions of the carrier density in samples SI, S2 and S3. Char- 
acteristic scattering times Ti and r* (in ps) are also shown for 
sample SI. 



voltage for the three samples. The bars indicate three 
regions of the gate voltage where p{T) was measured in 
all samples: V g = 3, 16 and 36 V with respect to the 
Dirac point. 

Figure^ shows p{T) in sample SI in the temperature 
range T= 5 - 200 K. The increase of the resistivity at high 
temperatures can be partially ascribed to the effect of 
acoustic phonon scattering on the classical conductivity 
Til 17 1, which is shown by the dashed line 18 1: 
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where D a =18 eV (as determined from the analysis of the 
classical conductivity [13]) is the deformation potential, 
p s = 7.6 x 10~ 7 kg m~ 2 is the density of graphene, v p h — 
2 x 10 4 m s _1 is the speed of sound, vf = 10 6 m s _1 is 
the Fermi velocity of carriers. We have subtracted the 
phonon contribution from the experimental dependence 
p(T). The analysis has been limited to the range T < 50 
K in order to rule out other types of phonons at higher 
temperatures [lil [l7| . The resulting quantum correction 
to the conductivity is shown for all regions in Fig. as 
Acr(T) = cr(T) - a {T ), where T is the lowest studied 
temperature, er(T) = [p{T) - p P h{T)]~ l . 

The separation of the EEI corrections from the WL 
contribution has been performed by two methods. For 
samples SI and S3 the low- field perpendicular magne- 
toresistance has been measured, in order to determine the 
characteristic times responsible for WL: the inelastic de- 
phasing time Tip(T), the elastic time of inter- valley scat- 
tering and the elastic time t* which describes intra- 
valley suppression of quantum interference (due to topo- 
logical defects and 'trigonal' warping of the energy spec- 
trum [l(|). (This analysis is done following the method 



described in |llLll2j.) These times are used to determine 
the WL correction 8a WL {T) 0, 
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which is then subtracted. In sample S2, the EEI correc- 
tion has been isolated by suppressing WL by a perpendic- 
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FIG. 2: (a) The resistivity as a function of temperature for 
sample SI for three regions of carrier density. The dashed line 
is the acoustic phonon contribution calculated using Eq.([2]) 
(the right-hand side axis), (b) The conductivity after the 
phonon contribution has been subtracted. Solid lines show 
the WL correction found from Eq.©. 



ular magnetic field which is still too small to affect the 
EEI correction [l[. Both methods lead to close results 
for the magnitude of the EEI correction in the studied 
samples. 

The solid line in Fig. shows the WL correction 
to the conductivity, A5a WL {T) = Sa WL {T) - <5cr WL (T ), 
found from the analysis of the magnetoresistance using 
the first method. One can see that the two types of 
quantum correction, WL and EEI, are of similar magni- 
tude. The solid lines show clearly that in regions I and II 
there is a transition from weak localisation, an increase 
of Acr(T), to antilocalisation, a decrease of Aa{T). (Ear- 
lier, such a transition was detected in the change of the 
sign of MR [I3|, with the transition temperatures of ~ 
10 K in region I and ~ 25 K in region II, which is in 
agreement with this experiment.) 

Figure shows the temperature dependence of the 
resistivity of sample S2 in the temperature range 0.25 - 
40 K. First, the phonon contribution, Eq.Q, shown by 
the dashed line is subtracted in regions II and III. The 
remaining quantum contribution to the conductivity is 
presented in Fig. [3b, for different magnetic fields. One 
can see that with increasing B there is a decrease in the 
slope of the temperature dependence until a saturation is 
reached. This is a signature that the WL correction has 
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been suppressed while the EEI correction is not affected 
by magnetic field. 

Indeed, the suppression of WL is expected at fields 
which are much larger than the so-called 'transport' field 
Btr — fc/2elp, where l v is the mean free path [ijj]. For 
sample S2 the values of B tr are 120, 70 and 45 mT for 
regions I, II and III, respectively, and therefore it is not 
surprising that WL appears to be suppressed at B = 1 
T, Fig.|3}D. On the other hand, the effect of the magnetic 
field on the EEI correction is due to the Zeeman splitting 
of the triplet 'channel' and is expected at higher fields, 
g*\iBB > ksT [2|, where g* is the Lande g-factor, and 
fiB is the Bohr magneton. For the g-factor in graphene 
~2 [20[ and temperatures above 1 K, this condition is 
satisfied at fields higher than 1 T. 

The extracted EEI correction is shown for samples SI 
and S2 in Fig. SI where we also add the result for sample 
S3 in region I. It is indeed logarithmic in temperature, 
Eq.([T]), with close values of A, A — 0.5 - 0.8, for all three 
regions of the carrier density in the studied samples. 

To interpret the obtained value of A, we note that 
the theory distinguishes between the contributions 
from different quantum states of two interacting elec- 
trons, commonly referred to as 'channels'. The coefficient 
A takes the form A = 1 + c (1 - ln(l + F CT )/F CT ), where 
Fq is the Fermi-liquid constant. While the first term in 
this relation represents the universal contribution of the 
'singlet' channel, the second (Hartree) term describes the 
contributions of c 'triplet' channels. For example, in a 
single-valley 2D system (such as in GaAs) the coefficient 
c = 3 due to identical contributions of three spin triplet 
states. (When the this degeneracy is lifted by magnetic 
field , two components become suppressed, resulting in 
c = 1.) 

In two-valley 2D systems (e.g., in Si-MOSFETs 
the situation is more complicated. In the absence of inter- 
valley scattering, the valley index v — ± is a good quan- 
tum number. In this case the overall number of channels 
is 16, due to four-fold spin degeneracy of two interacting 
electrons and an additional four-fold degeneracy due to 
the two valleys. This gives the prefactor c = 15. This 
result also holds if the inter-valley scattering is weak, 
ksT h/r.i, i.e. when the typical electron energy is 
larger than the characteristic rate of inter-valley scatter- 
ing. However, at low temperatures, fcsT -C h/ri, strong 
inter-valley scattering mixes the valleys and A takes the 
same form as in the single- valley case. 

Unlike in Si-MOSFETs, in graphene the valley dynam- 
ics is governed by two characteristic times, the inter- 
valley scattering time Ti and intra- valley dephasing time 
t* [lOj. In our experiments, the inter- valley scattering 
rate h/rt is of the order of 3 K, while the intra- valley de- 
phasing rate is above 20 K (see Table Q}. Thus the inter- 
mediate regime, H/ri < ksT < h/r^, becomes possible. 
In this case, the channels with two electrons from differ- 
ent valleys give no contribution. (This situation is similar 
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FIG. 3: (a) The resistivity as a function of temperature 
for sample S2, shown for three regions. The dashed line 
is the acoustic phonon contribution calculated using Eq.(f2|) 
(the right-hand side axis), (b) The conductivity Acr(T) = 
<t(T) — t(To) at different magnetic fields (the contribution of 
acoustic phonons has been subtracted). 



to the one which occurs for universal conductance fluc- 
tuations in graphene [2l|.) As there are two spin states 
per electron and two states for two electrons in the same 
valley, there are eight remaining channels, one of which 
is both spin and valley singlet, so that c = 7. Thus we 
arrive at the following expression for A in Eq.([T]): 



A{F°) = 1 + 7(1- ln(l + W) 



(4) 



We have confirmed this analysis by standard diagram- 
matic calculations, where we used a common assumption 
that all channels except for the singlet are described by 
the same Fermi-liquid parameter. 

Using Eq.(j4|) and the experimental values of A, Fig. 
we find the values of Fq to be between -0.08 and -0.13. It 
is interesting to note that the value of Fq found in GaAs 
and Si systems at r s ~ 1 is between -0.15 and -0.2 0,0,0]. 
To explain this relatively low value of Fq ss —0.1 found 
in our experiments, we note that in graphene this con- 
stant is suppressed due to the chirality of charge carriers, 
which prevents large-angle electron-electron scattering. 
Indeed, in the non-chiral 2DEG, the constant Fq can be 
found by averaging the electron-electron scattering am- 
plitude over all possible scattering angles (see, e.g., 0): 
Fq = — v{U(\p — p'D). Here U(q) is the Fourier com- 
ponent of the interaction potential, and v is the density 
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FIG. 4: The electron-electron interaction correction to the 
conductivity, a) The results for sample 1 (circles) are obtained 
by determining the WL correction using Eq.© (in region I 
the results for sample S3 are also displayed by triangles), b) 
For sample S2 the WL contribution is suppressed by magnetic 
field. Solid lines are fits to Eq. ([!]). 



of electron states per spin/valley. In a chiral system, the 
scattering amplitude for each electron is suppressed by 
the factor cos(0/2), where 8 is the scattering angle (see, 
e.g. [H), so that F °" = -u{U(\p - p'|) cos 2 ((9/2)). 

For a simple estimate away from the Dirac point, we 
use the Thomas-Fermi approximation for the interaction 
potential, U(q) ~ 2-Ke\j (q + k), with the effective charge, 
e\ = 2e 2 /(e+l) that includes suppression of the Coulomb 
interaction by Si02 substrate, e = 3.9. The screening 
parameter n = 4 • 2~Kve\ includes contributions from four 
degenerate single-electron states (in graphene the density 
of electron states is v — ep/2irvp). This gives 
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d0 



TT 2 sin § + 2a 
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where a = e^/fiVF ~ 0.88 is the dimensionless interaction 
constant (it is related to the parameter r s used in [3j-l8j as 
r s = \/2a). Evaluating this integral, we find Fq = —0.10, 
which is in agreement with our measurements. Note, 
that a similar calculation for a non-chiral electron liquid 
with two valleys gives a larger value Fq w —0.19 for 
the same value of a. Approximation ([5]) which neglects 
effects of strong interaction, such as the Fermi velocity 



and Z factor renormalisations 22| , is expected to be valid 
for a < 1, which is the case for graphene. Our result is 
in agreement with the value of Fq in [22| for the studied 
range of charge densities. (To compare Fq with Fq in 
[22}, one has to take into account that these quantities 
are related as Fq = 2Fq .) 

In summary, we show that electron-electron interaction 
plays an important role in the low-temperature conduc- 
tivity of carriers in graphene. Unexpectedly for the EEI 
correction, its magnitude is affected by the intra-valley 
decoherence rate due to elastic scattering. We find the 
value of the interaction parameter Fq in graphene, which 
is lower than in other 2D systems studied earlier. 
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